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Abstract. Aerosol Optical Depth (AOD) and Angstrom ex- 
ponent (AE) values derived with the MODIS retrieval al- 
gorithm over land (Collection 5) are compared with ground 
based sun photometer measurements at eleven sites spanning 
the globe. Although, in general, total AOD compares well 
at these sites (R 2 * values generally over 0.8), there are cases 
(from 2 to 67% of the measurements depending on the site) 
where MODIS clearly retrieves the wrong spectral depen- 
dence, and hence, an unrealistic AE value. Some of these 
poor AE retrievals are due to the aerosol signal being too 
small (total AOD < 0.3) but in other cases the AOD should 
have been high enough to derive accurate AE. However, in 
these cases, MODIS indicates AE values close to 0.6 and 
zero fine model weighting (FMW), i.e. dust model provides 
the best fitting to the MODIS observed reflectance. Yet, 
according to evidence from the collocated sun photometer 
measurements and backtrajectory analyses, there should be 
no dust present. This indicates that the assumptions about 
aerosol model and surface properties made by the MODIS 
algorithm may have been incorrect. Here we focus on prob- 
lems related to parameterization of the land-surface optical 
properties in the algorithm, in particular the relationship be- 
tween the surface reflectance at 660 and 2130nm. The re- 
trieval assumes that there is a linear equation that relates the 
reflectance in these two channels, with the value of the slope 
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(slope 66 o/ 2 i 3 o) determined, in part, by the infrared Normal- 
ized Difference Vegetation Index, (NDVIswir)- However, 
the assumed dependence of the slope on the NDVIswir is 
not supported by a MODIS based surface albedo climatol- 
ogy. The use of a modified relationship based on the albedo 
data improves the AE retrieval at the studied sites. The in- 
crease in the AE agreement fraction between MODIS and 
AERONET measurements is between 3 and 22 percentage 
units depending on the site. These results indicate that the 
surface reflectance assumptions, especially the slope 660 / 2 i 30 
in the MODIS algorithm is the major reason for the inaccu- 
rate AE values and the flawed aerosol model combining in 
the retrieval. However, at some of these sites, the new rela- 
tionship slightly reduces the correlation between the MODIS 
and AERONET AOD. This decrease indicates that the com- 
bination of the assumed surface and aerosol properties still 
do not match the actual properties under investigation. 


1 Introduction 

Aerosols play an important role in the Earth’s atmosphere. 

They affect, for example, climate, radiation budget and cloud 
processes. The direct effect of aerosols on climate involves 
both scattering and absorption of radiation, while the indi- 
rect effect of aerosols is related to their ability to modify the 

optical and physical properties and thus also the lifetimes 
of clouds. Aerosols are ubiquitous and therefore satellite 
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measurements can provide global aerosol information. How- 
ever, because of the huge variety of aerosol and surface 
conditions, passive satellite algorithms are underdetermined 
in their information content. Therefore, to retrieve realis- 
tic aerosol properties, the algorithms require reasonable as- 
sumptions and constraints on aerosol properties and the con- 
tributions of the underlying surface to the observed radiance 
field. The quality of the assumptions is critical, especially 
over land where the surface reflectance can dominate the 
reflectance measured at the top of the atmosphere. Due to 
the stronger signal from the surface, the error in the derived 
aerosol optical depth (AOD) is typically 10 times larger than 
the error in the retrieved surface reflectance (Kaufman et 
al., 1997). Consequently, developing satellite retrieval algo- 
rithms is hard. 

One of the most widely used passive satellite instruments 
in aerosol remote sensing is the Moderate Resolution Imag- 
ing Spectroradiometer (MODIS, Salomonson et al., 1989). 
MODIS instruments are aboard the Terra and Aqua satellites. 
Having a swath width of 2300 km. Terra MODIS and Aqua 
MODIS cover the Earth’s surface every 1 to 2 days. Begin- 
ning in early 2007, a consistent algorithm (Levy et al., 2007a) 
was used to process the entire MODIS data record, provid- 
ing the so-called Collection 5 (hereafter C5) data set. Levy et 
al. (2007a, and the on-line ATBD) have provided a detailed 
documentation of the algorithm and the changes from previ- 
ous versions. 

The C5 over land algorithm differs from prior versions 
in its retrieval process logic, as well as in its assump- 
tions about aerosol model optical properties and surface re- 
flectance. The algorithm makes a simultaneous inversion of 
MODIS -observed reflectance in three channels (0.47, 0.66 
and 2.13 pm), to retrieve total aerosol loading, proportion of 
fine/coarse aerosol particles, and surface reflectance. The 
products are three semi-independent parameters, the total 
AOD (or r) at 0.55 pm, the nondust or fine model weight- 
ing (LMW), and the surface reflectance at 2.13 qm. LMW is 
the fractional contribution of fine mode aerosol to the total 
AOD at 0.55 qm. 

C5 uses five aerosol models (Levy et al., 2007b), four 
of which are used in the standard algorithm: dust, weakly 
absorbing, moderately absorbing and strongly absorbing. 
Each of these four models has been determined by analy- 
sis of AERONET almucantar inversions (e.g. Dubovik et al., 
2002), and are comprised of two lognormal modes in their 
size distribution. The dust model is non-spherical and dom- 
inated by coarse mode, whereas the other three are spherical 
and dominated by the fine mode. The three fine models are 
different mainly in their single scattering albedo (SSA), such 
that the weakly, moderately and strongly absorbing models 
have approximate SSA at 550 nm of 0.95, 0.91 and 0.86, re- 
spectively. The dust model has a SSA of 0.95 (Levy et al., 
2010). The value of the FMW represents the “mixture” of 
one of the fine-dominated models (hereafter fine model) and 
the dust model, where the fine model is prescribed by geog- 


raphy and season. Thus, the FMW is not a weight between 
fine and coarse modes, but represents the proportion of non- 
dust aerosol model. A FMW value of zero represents the 
case where the dust model provides the best solution. The 
retrieved spectral AOD represents the optical contributions 
of the dust/non-dust mixture, from which the Angstrom ex- 
ponent (AE) can be derived. 

Validation with data from the Aerosol Robotic Network 
(AERONET, Holben et al., 1998) shows that retrieval of 
AOD has improved compared to C4 (Levy et al., 2007a; Mi 
et al., 2007; Papadimas et al., 2008; Jethva et al., 2007). 
Moreover, MODIS C5 measurements have been extensively 
compared with AERONET measurements globally (Levy et 
al., 2010; Levy et al., 2007a; Jethva et al., 2010; Jethva et 
al., 2007; Mi et al., 2007; Kaskaoutis et al., 2007; Papadi- 
mas et al., 2008; Misra et al, 2008; Drury et al., 2008; Re- 
nter et al., 2008, Mishchenko et al., 2010). Several studies, 
for example Jethva et al. (2010) and Oo et al. (2008) have 
concluded that the surface assumptions and aerosol models 
in the MODIS retrieval are not representative for all loca- 
tions. As for AE, there have been no comprehensive studies 
comparing MODIS and AERONET values of AE. Levy et 
al. (2007a) did not find high levels of correlation from their 
small sample of comparisons. Mishchenko et al. (2007), Liu 
and Mishchenko (2008) and Mishchenko et al. (2009) com- 
pared MODIS to retrievals from the Multi-angle Imaging 
SpectroRadiometer (MISR, Kahn et al., 2001) and found no 
correlation between the two instruments, further indicating 
that the information observed by passive satellite techniques 
over land is weak. Although Kahn et al. (2009) showed that 
some of the poor comparisons were resulting from data mis- 
interpretation, retrieval of global AE is still considered a dif- 
ficult problem. 

Why are the AE discrepancies so large, even when to- 
tal AOD is comparable? Why does MODIS retrieval report 
100% dust, even in regions, where no significant dust load- 
ing should be expected? This study addresses the question 
why the MODIS algorithm wrongly selects the dust model 
in such cases. Previous studies have suggested, and Levy et 
al. (2010) have stated, that the MODIS observations do not 
contain enough information to robustly separate dust from 
non-dust. Uncertainties in the surface reflectance, as well as 
in aerosol absorption are too large. However, from the collo- 
cations with sun photometer data, we can begin to quantita- 
tively evaluate these systematic errors. In some sites the mis- 
identification occurs because the surface reflectance assump- 
tions are consistently wrong, whereas at other sites, the error 
is dominated by the initial assumption of aerosol model type. 
In this paper, we study some of the sites where MODIS is 
consistently picking the wrong aerosol model, which in turn 
leads to errors in total AOD and especially in AE. We deter- 
mine the dominant cause of aerosol type mis-identification, 
and determine whether the retrieval can be tuned in a physi- 
cally meaningful way at these sites. 
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Table 1 . Sunphotometer sites used in the study. PFR sites were in Finland and AERONET data was used elsewhere. 


Site 

Country 

Lat., Lon. 

Elevation [m] 

PFR: 

Jokioinen 

Finland 

60.81° N, 23.50° E 

104 

Sodankyia 

Finland 

67.37° N, 26.63° E 

184 

AERONET: 
Alta Floresta 

Brazil 

9.87°S, 56.10°W 

277 

GSFC 

USA 

38.99° N, 76.84° W 

87 

Ispra 

Italy 

45.80° N, 8.63° E 

235 

Kanpur 

India 

26.51° N, 80.23° E 

123 

Mexico City 

Mexico 

19.33° N, 99.18° W 

2268 

Mongu 

Zambia 

15.25° S, 23.15° E 

1107 

Rome Tor Vergata 

Italy 

41.84° N, 12.65° E 

130 

Toravere 

Estonia 

58.26° N, 26.46° E 

70 

Xianghe 

China 

39.75° N, 1 16.96° E 

36 


2 Data and methodology 

In this study, AOD available from MODIS is compared with 
AOD measured with ground-based Precision Filter Radiome- 
ters (PFR, Kim et al., 2008) located in Finland. In addition, 
MODIS measurements are compared with Aerosol Robotic 
Network (AERONET, Holben et al., 1998) Cimel sun pho- 
tometer measurements. The sites are listed in Table 1. PFR 
measurements in Sodankyia and Jokioinen were made be- 
tween 2004 and 2008, and 2006 and 2008, respectively. For 
both sites, data from days with residual snow on the ground 
were removed from the analysis because snow would render 
the MODIS retrieval difficult. Cimel measurements at the 
AERONET sites were done between 2000 and 2008. 

2.1 MODIS data 

MODIS C5 data from both satellites are used in the compar- 
isons with the ground-based measurements. The expected 
error over land in the MODIS AOD (or r) is ±0.05±0.15r 
(Levy et al., 2010). We used level 2 AOD data which 
has a spatial resolution of 10 x 10 km 2 at nadir. Only 
AOD data with best quality (Quality flag = 3) were con- 
sidered in the analysis (Levy et al., 2010). For the com- 
parisons with sun photometer measurements we used Opti- 
cal_Depth_Land_And_Ocean product. In the more detailed 
analysis of the MODIS retrieval algorithm we used Cor- 
rected Optical Depth Land product. These data were down- 
loaded from LAADS Web (http://ladsweb.nascom.nasa.gov/ 
data/search. html). 

2.2 Sun photometer data 

Two different types of sun photometer instruments were used 
in this study. The Precision Filter Radiometer (PFR; Wehrli, 
2000; Kim et al., 2008) is a passive instrument which mea- 
sures direct solar irradiance in four narrow spectral bands 


(368, 412, 500, and 862 nm). The bandwidth of the instru- 
ment is 5 nm and the full field of view angle is 2.5 degrees. 
Derived products from the measurements are AOD for the 
four wavelengths and AE. The time resolution of the prod- 
ucts is 1 min. The absolute uncertainty in AOD of the PFR in- 
strument is between 0.01 and 0.02 (Carlund et al., 2003). For 
the comparison with MODIS, the AOD measured by the PFR 
at 500 nm was interpolated to 550 nm using the Angstrom 
power law, T550 = r5oo(550/500) _ “, where 7550 is the AOD 
at 550 nm, T500 is the AOD at 500 nm and a is the AE based 
on the linear fit taken over all four PFR wavelengths. Cloud 
screening of the PFR measurements is based on the method 
presented by Smirnov et al. (2000). 

The Aerosol Robotic Network (AERONET; Holben et al., 
1998) uses Cimel sun photometers which measure AOD at 
340, 380, 440, 500, 675, 870 and 1020 nm. Measurements 
are provided every 15 min during daytime. AERONET also 
provides the angular distribution of sky radiances at four 
wavelengths (440, 670, 870 and 1020 nm), and aerosol prop- 
erties such as aerosol size distribution, complex refractive in- 
dex, and single scattering albedo (SSA), once every hour in 
clear sky conditions. The direct-sun AOD from AERONET 
are accurate to within ±0.01 for wavelengths larger than 
400 nm and ±0.02 for shorter wavelengths (Holben et al., 
1998; Eck et al., 1999). Holben et al. (1998) have described 
AERONET measurements in more detail. For the compar- 
ison with MODIS, the AOD measured by AERONET at 
500 nm was interpolated to 550 nm by using the Angstrom 
power law. We used cloud screened and quality assured level 
2.0 data and level 2.0 inversion products in the analysis. 

We compared the AERONET (Cimel) and PFR AOD 
measurements carried out at Sodankyia between March and 
September, 2007. For 1790 coexistent measurements with 
less than 3 minute time difference, the comparison showed 
that the ground based instruments agree extremely well at 
all comparable wavelengths. At visual and UV wavelengths 
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Table 2. Statistics of the concurrent AERONET/PFR and MODIS AOD measurements. Subscript “all” refers to all coexistent measure- 
ments, “over” refers to measurements with MODIS AE over 1, and “under” refers to measurements with MODIS AE under 1. /z is the 
number of measurements. R 2 is the correlation coefficient, mean(gb-M) is the average difference between the ground-based and MODIS 
measurements, RMS diff is the Root-Mean-Square Difference, and Pixels is the percentage of pixels within the expected accuracy of MODIS. 
Coarse^ERONET an d Coarsejy[ODls are the percentages of coarse particle measurements at each site based on AERONET and MODIS data, 
respectively. 



Alta Floresta 

GSFC 

Ispra 

Jokioinen 

Kanpur 

Mexico City 

Mongu 

Rome 

Sodankyla 

Toravere 

Xianghe 

All measurements: 

"all 

244 

895 

652 

126 

340 

172 

1085 

785 

127 

252 

462 

*all 

0.97 

0.93 

0.90 

0.90 

0.80 

0.64 

0.87 

0.78 

0.80 

0.95 

0.94 

mean(gb-M) a n 

-0.02 

-0.02 

0.02 

0.01 

-0.06 

0.01 

0.01 

-0.02 

0.02 

0.02 

-0.02 

RMS diff a n 

0.20 

0.05 

0.06 

0.03 

0.14 

0.10 

0.08 

0.05 

0.03 

0.05 

0.12 

Pixels [%] 

67 

92 

86 

97 

79 

71 

85 

88 

92 

87 

80 

Coarse AERO net [%] 

14 

4 

6 

10 

49 

12 

4 

12 

0 

8 

30 

Coarse M QDIS [%] 

21 

40 

2 

37 

92 

67 

21 

70 

17 

16 

84 

Measurements with MODIS AE over 1 : 

w over 

193 

541 

640 

79 

26 

56 

862 

233 

105 

212 

76 

d2 

^over 

0.97 

0.96 

0.90 

0.90 

0.82 

0.68 

0.89 

0.77 

0.81 

0.96 

0.96 

mean(gb— M)over 

-0.03 

-0.02 

0.02 

0.01 

0.07 

0.03 

0.02 

-0.00 

0.02 

0.02 

0.02 

RMS diffover 

0.22 

0.04 

0.06 

0.03 

0.13 

0.08 

0.08 

0.04 

0.03 

0.05 

0.15 

Measurements with MODIS AE under 1 : 

n under 

51 

354 

12 

47 

314 

116 

223 

552 

22 

40 

386 

R^der 

0.64 

0.67 

0.72 

0.89 

0.82 

0.62 

0.56 

0.77 

0.80 

0.90 

0.93 

mean(gb-M) undcr 

0.05 

-0.03 

-0.00 

0.01 

-0.07 

0.01 

-0.01 

-0.02 

0.02 

0.00 

-0.03 

RMS diff undcr 

0.07 

0.06 

0.04 

0.03 

0.14 

0.10 

0.07 

0.06 

0.03 

0.04 

0.11 


(368, 412 and 500 nm) the correlation coefficients (R 2 ) val- 
ues were 0.99. For the longer wavelength band (862 ran) the 
correlation is almost as good, with R 2 = 0.97. For the all 
the wavelengths the bias between the instruments was less 
than 0.01. This and previous analyses (e.g. Kim et al., 2008) 
show that PFR and AERONET measured AOD are of similar 
quality. 

2.3 Collocation of MODIS and sun photometer data 

The ground-based measurements were used in the analysis 
only if the measurements were made within an one hour win- 
dow centered at the MODIS overpass. The averaging of the 
satellite and surface measurements was based on the method 
described by Ichoku et al. (2002). If MODIS had more than 5 
aerosol retrievals inside a grid of 5 x 5 pixels (50 x 50 km 2 ), 
centered at an AERONET/PFR station, these retrievals were 
averaged and compared to ground-based measurements. In 
addition, ground-based measurements were used in the anal- 
ysis only if there were more than 6 PFR or 3 AERONET 
measurements within the one hour window. This was done 
to ensure that one instrument did not measure in the begin- 
ning of the time window and the other one in the end. 

Comprehensive global comparisons between AERONET 
and MODIS are presented by Levy et al. (2010), here we 
compare MODIS and sun photometer data for the specific 
sites, where there are many cases where MODIS is not iden- 
tifying dust properly. Table 2 summarizes the differences be- 


tween MODIS and AERONET/PFR AOD measurements at 
the sites studied. All parameters are presented for all mea- 
surements, for the measurements with MODIS AE over 1, 
and for the measurements with MODIS AE under 1. The 
measurements were stratified by the MODIS AE in order to 
study how errors in AE are tied to errors in AOD. In addi- 
tion, Table 2 includes the percentage of coarse particle mea- 
surements at each site based on AERONET and MODIS data 
(amount of measurements with AE under 1). The AERONET 
values show that in Kanpur 49% of the coexistent measure- 
ments has coarse aerosols while in Xianghe the fraction is 
30%. The other sites do not have significant amounts of 
coarse aerosols. However, AE data from MODIS suggests 
that the coarse fraction is larger than 15% at all the sites, ex- 
cept in Ispra where the fraction is only 2%. In Kanpur the 
MODIS coarse fraction is as high as 92% while in Xianghe 
it is 84%. 

At Alta Floresta, GSFC, Ispra, Mexico City, Mongu, and 
Toravere the R 2 values are the lowest for the cases with AE 
under 1. At Jokioinen, Rome, Sodankyla, and Xianghe the 
values of R 2 for the measurements with AE under 1 are 
slightly lower or in the same range as in the other classes. At 
Kanpur the R 2 values are equal for the AE > 1 and AE < 1 
classes. The averaged difference between the instruments 
(ground based - MODIS, henceforth mean(gb-M)) and the 
Root-Mean-Square of the difference (RMS diff) are very 
similar for all categories. In addition. Table 2 shows the per- 
centage of the measurements within the expected MODIS 
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Kanpur 


Xianghe 


Kanpur 


Xianghe 






AERONET AOD 


AERONET AOD 


AERONET AOD 


AERONET AOD 


Fig. 1. Correlation of AERONET/PFR and MODIS measurements 
in Jokioinen, Mongu, Kanpur and Xianghe. Blue points indicate 
measurements with MODIS AE over 1, and red points measure- 
ments with MODIS AE smaller than 1. Statistics for the compar- 
isons are shown in Table 1 . 


Fig. 2. Correlation of AERONET/PFR and MODIS measurements 
in Jokioinen. Mongu, Kanpur and Xianghe. Blue points indicate 
measurements with AERONET AE over 1, and red points measure- 
ments with AERONET AE smaller than 1 . 


accuracy. At all sites the percentages are larger or in the 
same range as the results presented by Remer et al. (2008) 
and Jethva et al. (2007). Both reported that more than 70% of 
the measurements at 550 nm over land were within expected 
accuracy. 

Figure 1 shows examples of comparisons of the AODs 
from MODIS and ground-based instmments at 550 nm 
wavelength stratified by MODIS AE. For some sites (e.g. 
Jokioinen, Mongu), dust is known to be rare (see Table 2). 
However, it is common at those sites that MODIS retrieves 
AE < 1. The question is why? We found that: (1) Retrievals 
with AE < 1 are connected with AOD < 0.5. (2) MODIS re- 
trieves AE < 1 at all studied sites, thus the flawed dust model 
selection is not confined to specific aerosol types or surface 
types. For the sites where also AERONET indicated the pres- 
ence of coarse aerosol particles (e.g. Kanpur and Xianghe, 
Table 2) the MODIS AE is under 1 also for larger AOD val- 
ues and the R 2 values for the two size classes (under and over 
1) are equal or in the same range. 

Figure 2 shows the same data as Figure 1 but stratified 
by AERONET AE. Evidently, the number of measurements 
with AE < 1 is significantly lower than with the MODIS 
data at all sites. In addition, most of the measurements 
with AERONET AE < 1 are associated with very low AOD 
values. 

2.4 SSA studies 

We used single scattering albedo (SSA) data from 
AERONET level 2.0 inversion products to study whether ab- 
sorption could explain the differences between MODIS and 
AERONET AOD values. To avoid measurement uncertain- 
ties associated with low AOD values, only measurements 


with AOD at 550 nm over 0.3 (both MODIS and AERONET) 
were considered. The AERONET SSA value was thought 
to be representative for the MODIS retrieval if the time dif- 
ference between the measurements was less than 10 h. Fig- 
ure 3 shows two examples of comparisons of SSA spec- 
tra for MODIS a-priori aerosol models with SSA values 
from the AERONET inversion categorized by MODIS AE. 
At Toravere, the a-priori SSA spectra in the MODIS algo- 
rithm are in agreement with the AERONET data for both 
size classes. The classification of the coarse aerosol parti- 
cles coincides with a strong dust plume advected from the 
Sahara on the 22 March 2007 (based on HYSPLIT backtra- 
jectories (Draxler and Rolph, 2010; Rolph, 2010), thus for 
these cases MODIS selects the aerosol models correctly. On 
the other hand, the SSA spectra measured by AERONET in 
Mexico City are lower than those for the a-priori model, at 
each wavelength. Similar behavior was also detected at Is- 
pra and Mongu. On the basis of AERONET AE, none of 
these sites were influenced by coarse particles during these 
measurements. This can also be seen from the spectral be- 
havior of the AERONET SSA. At Toravere, the averaged 
SSA of the cases with MODIS AE < 1 looks like the dust 
model, and the average of the cases with MODIS AE > 1 
looks like the fine model. At Mexico City, on the other hand, 
the spectral dependence of the AERONET SSA does not look 
like either of the MODIS assumed models. However, both 
cases have the same wavelength dependence, indicating more 
absorbing fine mode aerosols than assumed in the model. 
The only difference between the types is that on average the 
measurements where MODIS detects AE < 1 are associated 
with smaller AERONET SSA values. This raises the ques- 
tion whether the absorption characteristics of the fine mode 
aerosols could affect the selection of the aerosol model com- 
bination in the MODIS aerosol retrieval. However, before we 


www. atmos-meas-tech. net/4/201/201 1/ 


Atmos. Meas. Tech., 4, 201-214, 2011 


206 


T. Mielonen et al.: MODIS aerosol type selection 


Toravere 



Fig. 3. MODIS aerosol model SSAs compared with AERONET SSA products at two sites: Toravere and Mexico City. SSA spectra for 
a-priori dust and fine models (moderately absorbing) used in MODIS retrieval are shown in red and blue, respectively. The magenta line 
represents the averaged SSA spectrum from AERONET measurements for observations whose MODIS AE is over 1 (fine model). The green 
color is an average for observations whose MODIS AE is smaller than 1 (dust model). At Toravere the dust measurements are from a single 
dust episode (22 March 2007). Only measurements with AOD(550nm) values larger than 0.3 (MODIS and AERONET) are considered. 


can answer this question we have to study how the aerosol 
model combination is affected by surface parameterizations, 
see below. 

2.5 Surface reflectance studies 

In the next step, we investigated how the C5 surface re- 
flectance estimates influence the retrieval of AE. We used in- 
dividual pixels and divided the data into two classes based on 
MODIS AE as before. Then we calculated the monthly mean 
surface reflectances at 660 nm from all measurements and 
normalized the surface reflectance values with these mean 
values to remove seasonal variations. This was done to re- 
move the effects of natural variation in surface reflectance 
properties due to seasonal cycle of vegetation. Figure 4 
presents the normalized surface reflectances for measure- 
ments with MODIS AE over 1, and under 1, as a function of 
AERONET AE. The surface reflectance data was binned into 
0.1 wide AE classes. This example from Mongu is also rep- 
resentative for other sites which had a sufficient number of 
measurements from both classes (Alta Floresta, Kanpur, Xi- 
anghe). At these sites, measurements with MODIS AE over 
1 had larger normalized surface reflectance values than the 
measurements with MODIS AE under 1. This was also seen 
in the 2.13 pm surface reflectance data. Surface reflectance 
should not depend on the AE, as the AERONET AE data 


shows in Fig. 4 where the normalized surface reflectance for 
the MODIS measurements with AE > 1 is almost constant 
as a function of AERONET AE. However, in the MODIS re- 
trieval, larger particles (dust model) are associated with lower 
surface reflectance values than small particles (fine model). 
Based on these results, it seems that both aerosol absorp- 
tion and surface parameterizations could have effects on the 
aerosol model combination in the MODIS retrieval. In or- 
der to separate their effects, we have to study the retrieval 
algorithm in more detail. 

3 Analysis of the MODIS retrieval algorithm 

To determine why MODIS tends to wrongly identify dust in 
certain scenes, we need to focus on the details of the MODIS 
algorithm. Specifically, we try to determine which algorithm 
assumptions are violated, and why the retrieval behaves as it 
does. 

The MODIS aerosol retrieval assumes that the measured 
reflectance ( p *) at the top of the atmosphere (TOA) is 
composed of atmospheric ’’path radiance” (p a (9 ,Oo,(p)) and 
transmission of surface reflectance through the atmosphere: 

P*(e,o 0 ,(p) 

= p a (9,9o,cP) + F d (Oo)T(0)p s (0,eo,(l>)/{l-sp s (0,8o,(l>)) ( 1 ) 
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AERONETAE 


where 

slope 660/213 o = slope^f™ + 0.002© - 0.27, (5) 

yi n *-660/2t30 = — 0.000250 + 0.033 (6) 

slope 470 / 6 6o = 0.49, and (7) 

yint 4 7o/66o = 0.005 (8) 


where in turn 


slope 


NDVIswiR _ 
660/2130 — 


1 0.48 ; NDVIswir < 0.25 
0.48 + 0.2(NDVI S wir - 0.25) ; 0.25 < NDVIswir < 0.75 
0.58 ; NDVIswir > 0.75 


(9) 


Fig. 4. MODIS surface reflectance at 660 nm, normalized with 
monthly mean values as a function of AERONET AE. All available 
measurements from Mongu were used. Measurements with MODIS 
AE over (blue) or under (red) than 1 are shown separately. Only 
measurements with AOD (550 nm) values larger than 0.3 (MODIS 
and AERONET) are considered. In addition, standard deviations of 
the surface reflectances inside the bins are shown with bars. 


Here 9 is the view zenith angle, Oq is the solar zenith angle 
and <!> is the azimuth angle of the scattered radiation from the 
solar beam. Fj(6o) is the normalized downward flux for zero 
surface reflectance, equivalent to the total downward trans- 
mission, T ( 0 ) is the upward total transmission into the satel- 
lite field of view, and .v is the atmospheric backscattering ratio 
(Kaufman et al., 1997). 

Since values of p a (9,9o,(p ), Fd(9o), T(9), and s are un- 
known properties of the atmosphere, and the surface re- 
flectance is also unknown, Eq. (1) cannot be solved directly. 
Instead, the MODIS retrieval algorithm is designed to solve 
Eq. (1), by referring to the lookup tables (LUT) contain- 
ing simulated values of p a (9,9o,<p ), Fd(9o), T(9 ), and s. 
Therefore, Eq. (1) can be inverted, solving for the surface 
reflectance, e.g.. 


p\e,e 0 ,(t>) 


p a (0,9o,</>) — p*(9,9p,<l>) 

^ (p a (9,9o,4>)-p*(0,9o,4))-F d (6 0 ) T (9) 


( 2 ) 


Starting with the 2.13 pm wavelength, the MODIS algorithm 
determines the surface reflectance that solves the equation, 
repeating for each of seven aerosol loadings (correspond- 
ing to AOD at 0.55 pm of 0.0, 0.25, 0.5, 1.0, 2.0, 3.0, and 
5.0), and for the fine and dust models, separately. For each 
loading/model scenario, the calculated surface reflectance at 
2.13 pm is transferred to 660 nm and 470 nm with the follow- 
ing Eqs. (Levy et al., 2007a): 


7*660 — 7*2130 ' slo P e 660/2130 + y int 660/2130 ( 3 ) 

and 


7*470 — 7*660 ' slo P e 470/660 + y int 470/660 ( 4 ) 


and 

© = cos _1 (— cos9qcos9 + sin0osin0cos0). (10) 

Determined by empirical fitting (Levy et al., 2007a) the 
slope 660 / 2 i 30 and y-int 660 / 2 i 30 > are coefficients to the lin- 
ear equation transferring surface reflectance from 2.13 pm 
to 660 nm; likewise similar coefficients transfer surface re- 
flectance from 660 to 470 nm. © is the scattering angle and 
NDVIswir is the Normalized Difference Vegetation Index at 
IR, defined as 

NDVIswir = ^' 24 ~ ^ 13 ■ (H) 

7*1.24 + 7*2.13 

where p\ u and p* 13 are the MODIS -measured reflectances 
at 1.24 pm and 2.13 pm. 

By these equations, the surface reflectance in the two vis- 
ible channels (470 and 660 nm) can be estimated from the 
reflectance at 2.13 pm. Equation 1 is then applied to estimate 
TOA reflectance (p*) in the two visible wavelengths. Note 
that the result is the estimated TOA reflectance in three chan- 
nels (470 nm, 660 nm and 2.1 pm), corresponding to seven 
aerosol loadings, for the fine and coarse models, separately. 
These estimated TOA values, can be compared with the ob- 
served TOA reflectance. By interpolating between entries of 
aerosol loading, the algorithm finds the exact loading which 
satisfies equality between the estimated and observed TOA 
reflectance at 470 nm. 

Now, the algorithm can mix the fine and coarse models. 
By iterating the model mixing (FMW) through steps of 0.1, 
while simultaneously recalculating surface reflectance that 
satisfies the surface reflectance parameterization, the algo- 
rithm arrives at different TOA estimates for different mix- 
tures. Equality is required at 470 nm and 2.1 pm, however 
there will be differences between estimated and observed val- 
ues at 660 nm. The particular loading/mixture that minimizes 
the 660 nm difference is considered as the best solution, thus 
retrieving both the AOD and FMW. Since spectral AOD is 
consistently represented by the LUT, it is easy to report AOD 
at 470 and 660 nm, and derive AE from these values with the 
Angstrom power law. 
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Having AE of 0.6, the dust model has a flatter spectral de- 
pendence than any of the fine models (AE= 1.6). Thus, for 
a given loading index (AOD at 550 nm), the coarse (dust) 
model contributes more to the TOA reflectance at 660 nm 
than does the fine model. Since the algorithm has to find 
the best match at 660 nm, this means that for a given value 
of surface reflectance at 660 nm, the algorithm must either 
adjust the total aerosol loading or the relative contribution 
of the coarse model. If the surface target is “brighter” than 
expected via the surface parameterization, MODIS may ob- 
serve larger TOA reflectance than modeled by Eq. ( 1). In this 
case, MODIS will compensate by picking a dust only aerosol 
mixture. As a result, MODIS will indicate dust aerosol even 
when other observations do not indicate dust in the column. 

3.1 MODIS standalone algorithm 

In order to trace the performance of the MODIS retrieval 
algorithm, we obtained a “standalone” version of the code. 
By standalone, we mean a version of the C5 process code 
that includes all the physics and mathematics of the oper- 
ational aerosol retrieval, without the subroutines related to 
cloud screening, pixel selection, or any of the other MODIS 
processing overhead (making HDF files, etc.). For inputs, the 
standalone code requires only values of the observation ge- 
ometry (angles, location, month) and measured spectral TOA 
reflectance. Its outputs include only AOD, AE and FMW, 
plus values of the 660 nm matching error. Therefore, the 
standalone version enables us to study the mechanics of the 
retrieval, and test how changes in assumptions can affect the 
retrieved parameters. Specifically, the standalone algorithm 
allows us to easily modify the surface reflectance parameter- 
ization and determine how changes affect retrievals of total 
AOD and AE. 

3.2 Evaluation of the surface assumptions 

The surface parameterization of Levy et al. (2007a), given 
by Eqs. (3-9), includes a dependence on the NDVIswir - 
(Karnieli et al., 2001). This vegetation index is similar to the 
standard Normalized Difference Vegetation Index (NDVI; 
red and NIR; Tucker, 1979) except the use of the two longer 
wavelengths (2.13 and 1.24 pm) helps to reduce aerosol con- 
tamination. Basically, the surface parameterization assumes 
that the slope660/2i30 increases as a function of surface- 
greenness. 

To evaluate the quality of the assumed surface properties 
in the algorithm we used a MODIS based surface albedo cli- 
matology created by Moody et al. (2005). The global cli- 
matology contains spatially complete white-sky albedo (bi- 
hemispherical reflectance) data from the years 2000-2004. 
Moody et al. have also made a five-year aggregate of high 
quality MOD43B3 data which was filled using an interpola- 
tion technique (accuracy of 3-8%; Moody et al., 2008). This 
aggregate is assumed to be representative of an average year. 


In addition, the data set contains surface albedo statistics (av- 
erage and standard deviation) sorted by ecosystems. These 
are available at various resolutions of each 16-day time pe- 
riod’s snow-free filled albedo maps for each separate year as 
well as for the five-year climatology. The statistics are avail- 
able at various resolutions between 0.5 x 0.5 and 90 x 360 
(degrees of latitude x degrees of longitude). For our com- 
parisons we selected the albedo statistics for the five-year cli- 
matology with the 60 x 120 resolution. This data provided us 
averaged albedo values at three wavelengths (0.66, 1.24 and 
2.13 pm) categorized by ecotypes at 9 different areas. Com- 
bined, these 9 areas cover the whole globe. With the 16-day 
time resolution we get 23 surface albedo values representing 
the yearly cycle for each ecotype in each of the 9 areas. The 
overall absolute accuracy of the MODIS surface albedo data 
is within 0.05 (Liu et al., 2009). 

Surface albedo is the ratio of the radiant flux reflected from 
a unit surface area into the whole hemisphere to the incident 
radiant flux of hemispherical angular extent (Schaepman- 
Strub et al., 2006), thus it is not view angle dependent. The 
reflectance in the MODIS retrieval, however, depends on the 
solar and viewing geometry. To minimize the discrepancies 
resulting from the comparison of slightly different quanti- 
ties we decided to compare ratios instead of absolute values. 
The selected ratios were the slope^j^j™ and NDVIswir- 

These parameters were selected because the slope^j^j^® 
is the most important parameter governing the wavelength 
dependence of the surface reflectance in the MODIS re- 
trieval and it depends on NDVIswir as shown in Eq. 9. The 
slope 660 / 2 i 30 also depends on the scattering angle (as can be 
seen from Eq. 5), however, the effect of the scattering angle 
is not as large as that of the NDVIswir- In addition to the 
slope66o/2i30- the retrieval algorithm uses a ratio of 470 nm 
and 660 nm surface reflectance (slope47o/66o)- However, this 
second slope does not have a strong effect on the retrieved 
fine model weighting because Rayleigh optical depth domi- 
nates the signal at the shorter wavelengths. 

3.3 Analysis with the standalone algorithm 

We carried out a number of sensitivity studies with the 
stand-alone version of the MODIS aerosol retrieval algo- 
rithm. We modified the surface reflectance parameters 
(slope 66 o/ 2 i 30 and slope 47 o/ 66 o) to check how they affect 
the aerosol model combination (i.e. FMW), AE, and overall 
AOD. The FMW parameter was most sensitive to the changes 
in the slope 660 / 2 i 30 values. Because the slope 660 / 2 i 30 de- 
pends mainly on the slope^^™ , we decided to concen- 
trate on it. Figure 5 shows a typical example from Mongu 
how the modification of the slope^J^^ affects the aerosol 
model combination. The black cross is the measured top of 
the atmosphere (TOA) reflectance at 660 nm while the ma- 
genta square is the result for FMW = 0.5 from the original 
retrieval. The red and blue markers are results (FMW = 0.5) 
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Fig. 5. Sensitivity of the retrieved TOA reflectance at 660 nm on 
changes in the slopej^^i^^ • The black cross is the observed TOA 
reflectance at 660 nm and the magenta square is the result from the 
original retrieval (FMW = 0.5). The red and the blue markers are re- 
sults from retrievals (FMW = 0.5) where the slope^^j™ is 10% 
smaller (0.9) and 10% larger (1.1), respectively. The errorbars show 
how much the combination of the dust and fine aerosol models af- 
fects the results. The top of the bar indicates unmixed dust model 
(FMW = 0.0) while the bottom of the bar is for the unmixed fine 
model (FMW =1.0). 


slope*0.9 i- 
measured reflectance 
slope'1.1 ^ 
original retrieval t 


from retrievals where the slope^)^™ ^ as t> een multiplied 
by 0.9 and 1.1, respectively. Modification of the slope does 
not affect the surface reflectance at 2.13 pm, whereas an in- 
crease of 10% in the slope increases surface reflectance at 
660 nm and 470 nm by 10%. The errorbars at 660 nm show 
how much combining of the fine and dust models affects the 
results. The top of the bar indicates FMW of zero (dust 
model) while the bottom of the bar is for FMW of one 
(fine model). In this case, the original retrieval resulted in 
the best fit with the unmixed dust model because the mea- 
sured TOA reflectance at 660 nm was larger than the calcu- 
lated one. When the slope^^Kio* was 10% larger than the 
original, the calculated TOA reflectance at 660 nm became 
larger and closer to the measured value, thus changing the 
selected aerosol model from dust to fine. On the other hand, 
if the slope^Y*^® was 10% smaller than the original, the 
calculated TOA reflectance at 660 nm moved further away 
from the measured value without affecting the aerosol type 
selection. 

As illustrated in Figure 5, the combination of aerosol 
models (and therefore AE) is sensitive to the value of the 
slope^^j^Q R , therefore it is extremely important that the 
slope value represents the surface under investigation very 
accurately. We made a sensitivity study which showed that 
in order to have less than 20% error in the retrieved AE value, 
the error in the slope^Y*^ 11 has to be less than 10% when 
AOD is between 0.2 and 0.4. For larger AODs the error in 
the AE is smaller but for smaller AODs it can be signifi- 
cantly higher. The error of 20% in AE means that AE of 


0.6 could vary between 0.48 and 0.72, and AE of 1.5 could 
vary between 1.2 and 1.8. This accuracy would still enable 
the separation between fine and coarse aerosol particles. 

To study the quality of the slope-NDVI relationship 
(slope^J^j™ in Eq. 9) in the algorithm, we compared 
it with the MODIS surface albedo climatology. Figure 6 
presents the dependence of the slope^j^ 1 ]™ on NDVIswir 
based on global MODIS surface albedo data. The different 
ecotypes are separated by colors and symbols. In the sur- 
face albedo climatology, the globe was divided into 9 areas 
(60 x 120 degrees each) and averaged surface albedo values 
for all ecotypes in these areas were calculated. Each ecotype 
from each of the 9 areas has 23 points representing the yearly 
cycle of the vegetation. The solid black line in each of these 
figures shows how the relationship of the parameters is taken 
into account in the MODIS aerosol retrieval. The albedo data 
in Figure 6 show that the relationship between NDVIswir 
and slope^^™Q R has a large variation between different 

ecotypes. However, for all ecotypes the slope^^*® de- 
creases as NDVIswir increases which is opposite to the re- 
lationship in the retrieval. This behavior was confirmed with 
linear regression models which all had negative slope values 
for each ecotype. In addition, the retrieval assumes smaller 
slope^j^Jj™ values than the surface albedo data shows. 

To test how the AE and AOD results are affected if we as- 
sume the slope-NDVI relationship to behave more like the 
albedo data suggest, we applied, as a first try, an inverse re- 
lationship that fits the data slightly better (shown in Fig. 6 
with dashed lines). As one can see, the inverse relationship 
fits well to number of ecotypes, e.g. savanna and shrub- 
land. However, it is not able to represent, for example, ur- 
ban, wetland and needleleaf ecotypes. In an attempt to eval- 
uate the MODIS AE quality we used AERONET AE data. 
The instruments were in agreement on the particle size if 
both AEs were either larger than 1.4 or smaller than 0.9. 
In other cases the instruments disagreed. The measurements 
with AEs between 0.9 and 1.4 were left out of the analy- 
sis to avoid the difficult comparison of mixed aerosols. In 
this study, we compared single measurements, whose time 
difference was less than 10 min. The number of usable mea- 
surement (Table 3) varies for the modified retrievals due to 
the “AE buffer zone” (from 0.9 to 1.4) and the exclusion 
of retrievals with AOD less than 0.3 or negative surface re- 
flectance at 2.13 pm. The comparison results are presented 
in Fig. 7 and in Table 3, which present the AE and AOD 
agreement parameters for 6 sites. The parameter AE agree- 
ment refers to the fraction of AE agreement while mean(gb- 
M) is the averaged AOD difference between the instruments 
and R 2 is the correlation coefficient for AOD. Retrievals 
done with the original slope-NDVI relationship as used in 
the C5 MODIS aerosol retrieval are called “Original” while 
the retrievals done with the inverse of that relationship are 
called “Inverse”. Both of these main retrievals have sub- 
retrievals where the fine aerosol model was assigned to be 
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+ Evergreen Needleleaf Forest 
+ Evergreen Broadleaf Forest 
Deciduous Needleleaf Forest 
+ Deciduous Broadleaf Forest 
Mixed Forest 
o Closed Shrubland 
o Open Shrubland 
Woody Savannas 
o Savannas 
Grasslands 

+ Permanent Wetlands 
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Urban and Built-Up 
+ Crop/Natural Vegetation Mosaic 


Fig. 6. Dependence of slope^J^**Q R on NDVI^wir for different ecotypes based on MODIS surface albedo data. Plot a) presents forest 
ecotypes (crosses), plot (b) savanna ecotypes (circles) and plot (c) other ecotypes (stars). The solid black line indicates the relationship of the 
parameters used in the MODIS retrieval (slope^J^™) while the dashed line is the inverse version of it. The inverse relationship appears 
to be a better fit for a number of ecotypes. 


either strongly absorbing (strong, SSA = 0.86) or weakly ab- 
sorbing (weak, SSA = 0.95) regardless of the fine model orig- 
inally assigned for the location. Table 3 and Fig. 7 basically 
show how the selection of the slope-NDVI relationship and 
fine aerosol model affect the AE and AOD agreement at each 
site. 

Figure 7 and Table 3 show that at all studied sites (Alta 
Floresta, Kanpur, Mexico City, Mongu, Rome, Xianghe) the 
modification of the slope-NDVI relationship improved the 
AE agreement. For example, in Alta Floresta the original 
AE agreement fraction is 0.83 and with the inverse relation- 
ship the AE agreement fraction increases to 0.92. When 
the inverse relationship is combined with the weakly absorb- 
ing model, the AE agreement increases to 0.94. Unfortu- 
nately, the AOD correspondence at the studied sites does 
not improve as clearly. The R 2 values between MODIS and 
AERONET AOD (R2 in Fig. 7) only improve at Alta Floresta 
when using the inverse slope-NDVI relationship. In Kan- 
pur and Mexico City the use of the inverse relationship im- 
proves the correlation but even better correlation is achieved 
by merely using a different fine model. However, the aver- 
aged AOD difference (mean(gb-M)) decreases at all sites, 
except in Mongu and Xianghe. In Mongu and Xianghe all 
three parameters give mixed results. The best AOD corre- 
spondence is achieved with the original slope-NDVI relation- 
ship while the best AE agreement is achieved with the inverse 
relationship. In Alta Floresta, AOD and AE correspondence 


can be improved at the same time by using the inverse of the 
slope and weakly absorbing model. The results for Rome 
are similar as for Alta Floresta, although the R 2 value for 
the Inverse(weak) case is not the best possible. In Kanpur 
and Mexico the two AOD correspondence parameters (R 2 
and mean(gb-M)) give mixed results. By using the inverted 
slope-NDVI relationship and strongly absorbing model in 
Kanpur we get the best values for AE agreement and mean 
AOD difference, but R 2 indicates the best AOD correspon- 
dence with the original slope-NDVI relationship and absorb- 
ing model. Jethva et al. (2010) also concluded that the 
MODIS retrievals over India could be improved by using a 
more absorbing fine model and increased surface reflectance 
in the visible channels. These results show that the absorp- 
tion of the fine model does not affect the AE agreement, how- 
ever, it has a significant effect on the AOD correspondence 
at some sites, like Alta Floresta and Xianghe. The results in 
Mexico City and Rome are poor in all the cases because they 
are urban sites with bright surfaces. In addition, Mexico City 
is an elevated site which also increases the uncertainty in the 
MODIS retrieval (Levy et al., 2010). Therefore, we stud- 
ied how much the slope^^S™Q R should be increased at these 
two sites in order to improve the AE and AOD retrievals. We 
did the retrievals several times and for each new run we in- 
creased the slope by 10%. Then we compared these new 
retrieval values with AERONET measurements and found 
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Table 3. Comparison of the AE and AOD agreement parameters from the original retrieval (Original), and the retrieval with the inverse 
NDVI-slope relationship (Inverse) at 6 sites. In addition, retrievals where aerosol models were assigned to be either strongly absorbing 
(strong, SSA 55 o nm = 0.86) or weakly absorbing (weak, SSA 55 o nm = 0.95) are presented for both classes. AE agreement refers to the fraction 
of AE agreement while mean(gb-M) is the averaged AOD difference between the instruments (AERONET-MODIS). R 2 is the correlation 
coefficient for AOD. These values are also presented in Fig. 6. The best values for each of the parameters are shown in bold. In addition, the 
number of measurements used in the comparisons are presented. 



Alta Floresta 

Kanpur 

Mexico City 

Mongu 

Rome 

Xianghe 

AE agreement: 
Original 

0.825 

0.903 

0.136 

0.804 

0.319 

0.787 

Original (strong) 

0.869 

0.910 

0.117 

0.804 

0.304 

0.804 

Original (weak) 

0.865 

0.917 

0.115 

0.803 

0.309 

0.792 

Inverse 

0.918 

0.940 

0.141 

0.835 

0.534 

0.875 

Inverse (strong) 

0.930 

0.945 

0.167 

0.844 

0.472 

0.872 

Inverse (weak) 

0.938 

0.929 

0.167 

0.823 

0.534 

0.837 

mean(gb-M): 

Original 

-0.032 

-0.128 

-0.192 

0.054 

-0.108 

0.006 

Original (strong) 

-0.352 

-0.125 

-0.201 

0.050 

-0.108 

-0.048 

Original (weak) 

-0.030 

-0.132 

-0.202 

0.054 

-0.109 

-0.007 

Inverse 

-0.006 

-0.057 

-0.145 

0.098 

-0.049 

0.125 

Inverse (strong) 

-0.309 

-0.051 

-0.131 

0.105 

-0.056 

0.018 

Inverse (weak) 

0.000 

-0.057 

-0.133 

0.104 

-0.049 

0.122 

R 2 : 

Original 

0.713 

0.795 

0.354 

0.638 

0.529 

0.765 

Original (strong) 

0.764 

0.802 

0.518 

0.630 

0.502 

0.694 

Original (weak) 

0.868 

0.799 

0.420 

0.625 

0.505 

0.933 

Inverse 

0.748 

0.788 

0.466 

0.517 

0.505 

0.715 

Inverse (strong) 

0.762 

0.800 

0.454 

0.520 

0.498 

0.647 

Inverse (weak) 

0.873 

0.783 

0.378 

0.514 

0.505 

0.768 

Number of measurements : 
Original 

97 

217 

81 

214 

116 

47 

Original (strong) 

111 

221 

77 

230 

112 

56 

Original (weak) 

122 

217 

78 

233 

110 

53 

Inverse 

110 

184 

64 

170 

73 

40 

Inverse (strong) 

129 

183 

66 

167 

72 

47 

Inverse (weak) 

129 

183 

66 

164 

73 

43 


that in Mexico City the slope^^™ should be as much as 
1 .6 times larger than the original to introduce significant im- 
provements into the AE (AE agreement fraction 0.80) and 
AOD (R 2 = 0.56 and mean(gb-M) = —0.03) agreements. In 
Rome, the best AE (AE agreement fraction 0.95) and AOD 
agreement (R 2 = 0.57 and mean(gb-M) = —0.01) is achieved 
with 30% larger slope Similar results were pre- 
sented by de Almeida et al. (2007) who analyzed the ef- 
fect of the surface reflectance ratio between the 660 nm and 
2.13 pm on AOD retrievals over Mexico City. They found 
that by increasing the 660/2130 ratio from approximately 
0.55 to 0.73 the agreement between the MODIS and sun pho- 
tometer AOD improved significantly. 


4 Discussion and Conclusions 

We compared Collection 5 MODIS data to PFR and Cimel 
(AERONET) measurements at eleven locations in Europe, 
Asia, Africa North America and South America. In north- 
ern Europe, the AOD values were quite low ( <0.3), espe- 
cially in northern Finland. AOD values from MODIS and 
PFR/AERONET in Kanpur, Rome, and Sodankyla agree rea- 
sonably well (R 2 between 0.78 and 0.80) while AERONET 
and MODIS measurements in Alta Floresta, GSFC, Ispra, 
Jokioinen, Mongu, Toravere, and Xianghe agree very well 
(R 2 over 0.87). In Mexico City the R 2 value was only 0.64. 

The MODIS retrieval algorithm over land often selects the 
dust aerosol model for locations where dust is unlikely to be 
found. In addition, there seems to be very little combining 
of the fine and dust models, meaning that the FMW is typ- 
ically 1 or 0. The selection of the dust model reduced the 
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Fig. 7. Comparison of the AE and AOD agreement parameters from the original retrieval (Original), and the retrieval with the inverse 
NDVI-slope relationship (Inverse) at 6 sites. In addition, retrievals where aerosol models were assigned to be either a strongly absorbing 
(strong, SSA 55 o nm = 0.86) or weakly absorbing (weak, SSA 55 o nm = 0.95) are presented for both classes. Colors of the bars refer to these 
different retrievals. AE agreement (top panel) refers to the fraction of AE agreement while mean(gb-M) (middle panel) is the averaged AOD 
difference between the instruments and R2 (bottom panel) is the correlation coefficient for AOD. These values are also presented in Table 3. 


correlation between MODIS and ground based AOD at the 
sites where no significant amount of dust is expected. This 
happened especially when AOD values were small. Accord- 
ing to Levy et al. (2010) MODIS does not provide quanti- 
tative information about aerosol size (AE) over land. How- 
ever, it is an important parameter if one tries to estimate the 
anthropogenic component of aerosols. 

Our findings are in agreement with the recent studies pre- 
sented by Levy et al. (2010) and Jethva et al. (2010). In 
contrast to their studies, we also analyzed single MODIS re- 
trievals and therefore were able to study the combining of 
the fine and dust models in more detail. When employing the 
spatio-technique of Ichoku et al. (2002), it is harder to see 
what happens in the combining of the aerosol models due to 
the averaging of at least 5 measurements. Single measure- 
ments, on the other hand, show clearly that the retrieval al- 
gorithm usually selects either fine model weighting of 1 or 0. 
Because AE depends only on the combination of the models, 
the AE retrieval fails if the combining fails. Our study sug- 
gests that the surface reflectance assumptions, especially the 
slope^^j^ in the MODIS algorithm is the major reason 
for the inaccurate AE values and the flawed aerosol model 
combining in the algorithm. Comparison of the surface pa- 
rameters shows that the slope-NDVI relationship used in the 
MODIS retrieval algorithm is not supported by the MODIS 
albedo climatology: 


1 . the albedo data show that the slope-NDVI relationships 
depend strongly on ecotypes. The possibility to take 
this into account in future Collections should be studied 
further. 

2. the slope^y*j^ R values based on the climatology were 
larger than the ones in the algorithm. 

3. the albedo data also show that the slope^^^i^o* de- 
creases as the NDVI increases for all ecotypes which 
is opposite to the relationship in the algorithm. 

4. changing the slope-NDVI relationship in the algorithm 
to agree better with the albedo data (inverse relation- 
ship) increased the correspondence of MODIS and 
AERONET AE at the studied sites. AOD correlation 
did not improve as clearly, but the mean difference de- 
creased at all sites, except in Mongu and Xianghe. 

5. the uncertainty of the fine model’s SSA does not have a 
significant effect on the AE agreement between MODIS 
and AERONET, however, it affects the AOD correspon- 
dence. 

Although, the modification of the slope^^!™ i m P rove d 
the AE retrieval at all studied sites and AOD retrieval at some 
sites, the effect on global scale is not yet known. This will be 
studied further in the future. 
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